from numpy import *
from scipy import *
def:
    dy= zeros(2,float)
    dy[0] = y[1]
    dy[1] = 1000*(1 - y[1]^2)*y[2] - y[1]
    return dy

